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C/3 . Abstract 

. ,-H . We take the Gauss' linking integral of two curves as a starting point to 

£^ ' discuss the connection between the equation of continuity and the inho- 

mogeneous Maxwell equations. Gauss' formula has been discussed before, 
as being derivable from the line integral of a magnetic field generated by a 
steady current flowing through a loop. We argue that a purely geometrical 
result - such as Gauss' formula - cannot be claimed to be derivable from a 
law of Nature, i.e., from one of Maxwell's equations, which is the departing 
point for the calculation of the magnetic field. We thus discuss anew the 
derivation of Gauss' formula, this time resting on Helmholtz's theorem for 
vector fields. Such a derivation, in turn, serves to shed light into the con- 
nection existing between a conservation law like charge conservation and 

^\ . the Maxwell equations. The key role played by the constitutive equations 

^^ ' in the construction of Maxwell's electromagnetism is briefly discussed, as 

f^ I well. 

•^ ; 1 Introduction 

r~| I Dated January 22, 1833, there is a note by Carl Friedrich Gauss that appeared 

P !■ in Volume V of his complete works, a note that refers to a formula which is 

nowadays known as the Gauss' linking integral. It represents the very first 
result of what has grown to be in the following centuries a new mathematical 
K> I theory: knot theory, a branch of topology. Gauss presented his result, a formula 

S , for calculating the linking number of two curves, without giving any proof of 

it. In modern notation, what Gauss asserted was the following. Consider two 
non-intersecting curves, Ci : x(s) and C2 ■ y{u), (see Fig. 1), and define the 
linking integral of the two curves, LkiCi, C'2), to be given by 
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Figure 1: Two non-intersecting curves Ci : x(s) and C2 : y{u), whose linking 
number is equal to one, in this case. 



where the prime denotes derivation with respect to the curve parameter. 

According to Gauss, when Ci and C2 are closed curves, Lk{Ci,C2) equals 
the number n of times that one curve winds around the other: Lk(Ci,C2) = n 
(amazingly, it is the factor l/47r what makes Lk{Ci, C2) an mief^er-valued func- 
tional of the loops). Furthermore, Lk(Ci, C2) is invariant under re-parametrization 
and under continuous deformation of the curves. 

In the course of time. Gauss' linking integral has been generalized within 
the context of differential geometry [1]. It has been defined for any two oriented 
submanifolds in an Euclidean n-space E", whose respective dimensions, r and 
s, have to satisfy r + s + 1 = n. The special case of two curves in an Euclidean 
3-dimensional space corresponds to r = s = 1 and 71 = 3. The generalization is 
based on a concept called the degree of a mapping. In order to introduce such 
a concept, one has to rely on several tools from modern differential geometry, 
and it is based on these tools that a formal proof of Gauss' assertion has been 
given. However, a formal demonstration of this sort seems to be completely un- 
illuminating. Surely, Gauss' steps leading to equation (1) must have followed a 
very different path. Unfortunately, it seems difficult to find in current literature 
an elementary derivation of equation (1). It should be mentioned that the 
editor of Gauss' Works conjectured that the method used by Gauss in arriving 
at equation (1) relied on calculating the line integral of the magnetic field due 
to a steady current. Such a conjecture seems to be substantiated by the fact 
that Gauss' linking integral is presented in the volume of his Works dedicated 
to mathematical physics, in the section of electromagnetism. However, Gauss' 
interest in arriving at a formula for the linking number seems to stem from his 
duties as an astronomer: he wanted to know which regions of the celestial sphere 
had to be observed over a period of time, in order to register the passing of an 
asteroid or other astronomical objects [2]. The Earth's orbit defines one curve, 
whereas the orbit of the celestial object defines a second one. The linking integral 
of these two curves can be shown to determine whether the celestial object can 



be observed or not. Within this context, a magnetic field does not play any role, 
of course, and a derivation of equation (1) from magnetostatics seems unlikely. 
Whether Gauss relied or not on magnetostatics for deriving his result will remain 
an open question probably forever. Anyhow, it has been recently published an 
illustrative and instructive article [3] in which the connection to magnetostatics 
is presented as an adequate starting point for deriving Gauss' assertion. 

Let us see in detail how the connection between magnetostatics and equation 
(1) arises. Consider the line integral ^B(y) • dy of a magnetic field around a 
closed curve, C2 : y{u). The magnetic field is assumed to be produced by a 
steady current / flowing through another loop, Ci : x(s). The field B can be 
obtained by using the well-known Biot-Savart formula [6], 

B(y)^^/«xxJl^. (2) 

^^ Jci |y-x| 

It leads to the following expression for the line integral to be calculated: 

lj,oI dxx (y-x) 



= ^ / / (y-x) ■ (rfyxdx) 
47r Jc, M |y - x|" 
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where the last step follows by applying the cyclic property of the triple product. 
Alternatively, we can calculate f„ B(y) • dy by using Ampere's law: 



/ B(y) 



dy ^ / V X B ■ dS = /io / J ■ dS = iiohnd = l^onl. (4) 
Js Js 



C2 Js Js 

Here, lend means the total amount of current crossing the surface S. Comparing 
equations (3) and (4) we obtain Gauss' result, equation (1) with Lk{Ci,C2) = 
n. We stress, however, that equation (1) involves nothing but geometrical (or 
topological) properties of the curves, whereas equation (3) and Ampere's law 
do express some intrinsic property of the magnetic field. 

What lies at the origin of the laws of Biot-Savart and Ampere are Maxwell's 
equations. These equations express a law of Nature. In principle, we should 
not take them as a starting point, in order to demonstrate a purely geometrical 
property of two curves. It is certainly appealing to consider that given a current 
in a loop, it produces a magnetic field which winds around the loop. But - let us 
stress it again - this follows as a consequence of a law of Nature, which rules how 
a magnetic field can be produced by the flow of electric charges. At most, the 
picture of a magnetic fleld winding around a circuit could have inspired Gauss 
to arrive at his assertion. At proving it, however, he must have followed quite 
another path, in accordance with the requirements of utmost rigor that he put on 



any mathematical result -'^. To some physicists, questions of mathematical rigor 
might appear superfluous or not worth to deal with, particularly in cases where 
one can arrive straightforwardly at a formula, by starting from well established 
natural laws. To others, however, it results obvious that a mathematical truth 
cannot follow from a law of Nature. This, for the very simple reason that 
what we call a natural law might change or become obsolete in view of new 
facts, whereas a mathematical result, once it has been proved, will remain true 
forever. 

A parable might be useful to illustrate our concern. Let us assume that the 
following series of events takes place: 

1. In 1785 Coulomb reports his findings about a law of Nature, which now 
bears his name: Fc ^ <Zi<Z2/^^- 

2. In 1833 Gauss reports a formula involving the sum of the angles of a 
spherical triangle. However, he does not include a proof of his formula. 

3. During the 20th century mathematicians develop differential geometry 
and by applying this theory they prove Gauss' formula, establishing it as a 
mathematical truth. 

4. In 1998 a physics journal publishes a paper in which it is claimed that 
Gauss' formula can be proved by simple means. Starting point of the proof is 
Coulomb's law. 

5. In 2020 very accurate experiments show beyond any reasonable doubt 
that Coulomb's inverse-square law is not exact: It should be replaced by a 
power law in which the exponent of r is —2.000000534(13). 

We are obviously entitled to ask: Was the proof published in the physics 
journal flawless? At first sight and in view of point 5, it was not. 

Now, let us put aside point 5 - or assume it never happens - but ask anew: 
Can the proof published in the physics journal be correct? The mere possibility 
that point 5 could happen makes the proof based on Coulomb's law at least 
suspicious of being flawed. In any case, it raises interesting questions: If the 
proof is correct, then it should be based only apparently on Coulomb's law. 
Perhaps it was some geometrical property lurking behind the inverse-square 
law what was actually used in the proof. Or perhaps there exists some deep 
connection between Coulomb's law and spherical symmetry that precludes point 
5 from ever occurring. Bringing such a connection into light would certainly 
help saving time and efforts dedicated to test the validity of Coulomb's law. For 
instance, sophisticated experiments are continually planned to test natural laws 
like the gravitational inverse-square law [5]. To be sure, as physicists we are 
ready to accept the verdict of an experiment. But we must be able to properly 
interpret what the verdict means and what it does not^. 



^As just two examples of this, let us remind Gauss' unpublished work on non-Euclidean 
geometry, that he developed 30 years before Bolyai and Lobachevsky, as well as his also 
unpublished work on special complex functions. These works remained unpublished because 
they did not satisfy the requirement of utmost rigor, that Gauss so fervently demanded from 
his and other's work. 

^As a matter of fact, a phenomenon like e.g. the Lamb shift can be interpreted as an indirect 
demonstration that Coulomb's law does not hold true at very small distances. Alternatively, 



One could still argue that we should distinguish between a physical law 
and the mathematical statement in terms of which we express that physical 
law. The physical law comes in when we assert that the real world behaves as 
described by the mathematical statement. From this mathematical statement 
we may deduce other such statements by a series of if-then relationships. In 
the case of Coulomb's law, such a relationship would tell that if the field of a 
point charge goes as 1/r^ then some other results follow. If some experiment 
then shows that the real world does not behave according to the 1/r^ law, it 
does not mean that the if-then relationships turn to be false. This is true, 
of course. However, without going into philosophical questions concerning the 
foundations of mathematics, we may see it as embracing different theories that 
can be constructed independently from one another, as free creations of the 
human mind. Each mathematical theory consists of a set of axioms, definitions, 
and theorems. It is clear that we cannot use a statement A of one theory - be 
it an axiom or a theorem - to prove a statement B of another theory, unless 
A is shared by both theories. Let one of these theories be the mathematical 
idealization of the real world that we may simply call, for the argument's sake, 
electrodynamics. Let the other theory be pure geometry. We could hardly say 
that the 1/r^ statement of Coulomb's law, or more generally, that Maxwell's 
equations fit into the body of pure geometry. They are neither axioms nor 
theorems of pure geometry and therefore cannot be used as starting points 
to prove purely geometrical statements, like those involving the angles of a 
spherical triangle or the properties of two curves. 

Let us write Coulomb's law in vectorial form, Fc ^ 9i(Z2(i"i — 1*2)/ |i"i — i"2| , 
and so we have the integrand of Biot-Savart's law showing up. To put Coulomb's 
law to a test is as justified as to put Biot-Savart's law to a test. Coulomb's law 
can be traced back to one of Maxwell's equations: V ■ E = p/eq. Biot-Savart's 
law can be traced back to V x B =/ioJ, a special, steady-state, case of one of 
the Maxwell's equations. Change the spherical triangle in the parable above 
by Gauss' linking integral and the foregoing remarks apply to the real case. 
This real case offers us the opportunity to analyze the content of Maxwell's 
equations from a somewhat new perspective. In the following, we will undertake 
an analysis of the Biot-Savart law and the Maxwell's equations, on which the 
proof of Gauss' linking integral that appeared in the real physics journal [3] was 
based. This will lead us, at the end, to find the connection between the principle 
of charge conservation and Maxwell's equations. 

2 Analysis of the Biot-Savart law 



The argument which allegedly proves Gauss' result was based on equations 
(3) and (4). Let us see the extent to which such a proof effectively depends on 
electromagnetism. Assuming that the derivation is correct, it should be possible 
to reformulate it in terms of purely geometrical concepts. In order to do this, 

it can be interpreted as showing the Hrnit of vahdity of a theory based on nonquantized fields. 



we first try to see the essence of the arguments: equation (3) derives from 
the expression that gives the magnetic field at point y produced by a steady 
current-density J(x): 

B(y) = ^ / J(x) X -y^£x. (5) 

' ^^Jv ^ ^ |y-x|^ ^ ' 

Indeed, by taking J as appropriate for the current in a loop, the volume 
integral in equation (5) reduces to a line integral {J{'x.)(Px -^ Idx): 

ic |y-x|' 



B(y) = ^ /„rfxx^:^-^, (6) 



which is equation (2), the expression commonly referred to as the Biot-Savart 
law [6]. 

Equation (5), in turn, follows as a solution of one of the Maxwell's equations 
(for steady fields): 

V X B =/ioJ, (7) 

together with V • B = 0, although this last equation is not essential for our 
present purposes'^. 

Let us recapitulate the essential steps leading to Gauss' result. We have 
as a starting point Maxwell's equation for a stationary current, equation (7). 
From it, we derive Ampere's law, §„ B(y) • dy =fio J„J ■ dS, as well as the 
Biot-Savart law, equation (6), which leads to equation (3). In Ampere's law we 
must put JgJ ■ dS =Ienci- Here, lend means the total current going through 
the loop C2 that encloses the surface S. In Biot-Savart's law, the current / 
refers instead to the current in the wire. When this wire winds n times around 
C2, we have lend ~ nl and Gauss' result follows. 

It seems therefore that the essential points in the foregoing proof of Gauss' 
result are: 1) Maxwell's equation (7), and 2) the concept of current as a flow of 
charged particles. Our question is: Is it possible to extract from these two points 
some purely geometrical properties on which to base an alternative proof? 

Let us analyze more closely point 2). After all, the current might be merely 
playing the role of a "counting device" , and so it might be an exchangeable unit 
within the foregoing reasoning. Alas, such a counting property should not be 
intimately tied to the physical nature of the current. The factor / appearing in 
equation (6) corresponds to the current flowing through a thin wire, the loop Ci 
in our case. The replacement 3{x)d^x -^ Idx leading from equation (5) to (6) 
seems to be crucial. It follows from the way we relate our concepts of current- 
density and current. Indeed, such a replacement comes from considering that 
the volume element d^x appearing in equation (5) can be written as d'^x = Sodx, 

•^Notc that equation (5) automatically satisfies V • B = 0. We remark that equation (5) is 
a special solution of (7), whereby the general solution can be obtained by adding to (5) a term 
Vf7, with U an arbitrary scalar function. For the purpose of demonstrating Gauss' result we 
could have used the general solution, because §^ V(7 ■ dy = 0, so that it is not essential to 
invoke V • B = 0. 



where So is the cross-sectional area of the wire and dx a distance-element along 
it (the direction of dx being given through a unit vector n, parallel to J). Thus, 
J(x)d^x = 3{x)Sodx = |J(x)| Sondx — Idx, whereby the last step expresses the 
connection between current-density and current: |J(x)| Sq — I. 

Now, when going from equation (5) to (6), the current / refers to the charge 
going through each piece dx of the wire. On the other hand, when we use 
equation (4) we must replace the right-hand side, Jg 3 ■ dS, by the total current 
lend crossing the surface S which is enclosed by C'2 ■ When the current-carrying 
loop crosses C2 a number n of times, we must put lend = /g J • dS = nl 
in Ampere's law. As we said before, replacing the left-hand side of equation 
(3) by iJ-oTT'I Eind cancelling the common factor iiqI, we obtain Gauss' result, 
Lk{Ci,C'2) = n. Thus, both the equation lend = nl and the relationship 
between current-density and current, as expressed by 3{x)d'^x = /dx, seem to 
be intimately tied to our physical concept of current as being a flow of particles. 
However, the cancelation of / in the last step of the proof seems to indicate that 
the current was not so essential. We will see in the following that we can in fact 
dispose of the concept of a current, and that we can replace Maxwell's equation 
(7) by another one which expresses a purely geometrical property. 

3 Derivation of the Gauss Linking Integral 

Let us now start by considering the closed curve Ci : x(s). This curve can be 
embedded in a family of curves, under quite general conditions. To define a 
family of curves amounts the same as to define a vector field v(r). Indeed, the 
curves pertaining to the family of curves can be seen as integral curves of the 
vector field, i.e., as curves whose tangent vectors coincide with v(r). They fulfill 
the following equation: 

We will be finally interested in only one member of the family of curves, 
namely on Ci. In case that all the curves of the family are closed curves, as we 
will assume henceforth, we have^ 

V • v(r) = 0. (9) 

Let us remark in passing that equation (9) holds true also for the velocity 
field of an incompressible fluid of constant and uniform density p. In absence of 
sources and sinks, such a velocity field satisfies the continuity equation 

-^ + V ■ pv = pV ■ V = 0. (10) 
at 



*We recall that the coordinate-free definition of the divergence of a vector field F(x) is 
given by V ■ F(x) = lim^'— .0 T7 /g -^ ' nda. When all integral curves of F(x) are closed, the 
surface integral in the definition of the divergence gives zero. 



Although in the present case equation (9) follows because the integral curves 
of V are closed, its relationship with the continuity equation will be important 
afterwards. Now, for v satisfying equation (9) and under appropriate conditions 
(see Helmholtz's theorem below), there exists a field W, whose curl is v: 

VxW = v. (11) 

This last equation is analogous to equation (7). Wc have so introduced the 
vector fields v and W, playing the roles of J and fJ.^^'B, respectively. From 
equation (11) we can derive the analogous to Ampere's and Biot-Savart' laws: 
equation (11) can be solved for W explicitly, by writing W = V x Z, with 

From this Z, by taking its curl, we obtain a solution of equation (11) given 

by 

W(y) = -^ / v(r) X ^^^d'r. (13) 

47r Jv \y - r\ 

It corresponds to equation (5) for the magnetic field. In order to derive Gauss' 
result from here, by following similar steps as before, we would need to reduce 
the volume integral in equation (13) to a line integral. To this end, we invoke 
the Dirac delta function to single out the contribution to W(y) which comes 
from a single integral curve x(s) of v(r). This is equivalent to restricting the 
field v(r) in equation (13) to be different from zero only along the curve x(s)^: 

= / ^4^,5-^(r - x(s))ds. (14) 

J as 

Introducing this expression for v(r) into equation (13) we obtain 



±1 U ,,!Mf), y-x(^) 

'^T^JcAJc, ds |y-x(s)| 



ds^^^ X , ' ";^/,3 . dy. (15) 



We have, therefore, that 



■''As in the commonly employed definition of a current density j(x, t) due to a single charge 
moving along the path r(t), i.e., j(x,t) = ev(t)5^(x — r(t)), an expression like equation (14) 
makes sense only under the assumption that subsequent integrations over its arguments have 
to be carried out. Equation (14) is similar to the covariant form of the 4-vector current 
density j°'{x), given in terms of the path r{T) that is traced back by a charge in space-time: 
j"(z) = ec/drj)°(r)5('*)(x-r(T)). 



This last expression can in turn be changed as follows, by applying the cyclic 
property of the triple product: 



£wiy).dy^l- 



1 

1 


i 
i 


/ 

/ 

JC2 


(x-y) 

X 

(y(«) - 


■ (dx X 

-yf 

x(s))- 


dy) 
iy'iu) X 


x'(s)) 


47r 




|y(") - 


-xGs)|^ 





'C2 

-duds. (17) 

We have so obtained the right-hand side of Gauss' formula. In this case we 
have made no reference to a current. As a second step, let us turn our attention 
to the left-hand side of equation (17). We have, by applying Stokes theorem 
and considering equation (11), that 

f W(y) ■dy= / (V X W) • iidS = / v • ndS, (18) 

Jc2 Js Js 

S being the area enclosed by C2 ■ We have to prove that the integral /„ v • ndS, 
with V given by equation (14), readily equals n, the number of times that one 
curve winds around the other, i.e., we have to demonstrate J„v ■ ndS = n. If 
we succeed, we have the aforementioned "counting device" at our disposal, and 
with it the second piece that is needed for the proof. Because the proof that 
J„ V • xidS = n is somewhat involved, we relegate it to an Appendix. 

Alternatively, we may introduce a small change in what we have done so far, 
in order to give a proof of Gauss' result that resembles closely the one based 
on magnetostatics. However, it is based on the concept of the strength of a 
tube of curves, a concept whose properties might be difficult to substantiate by 
purely geometrical means. It makes possible to show that J„ v • xidS works as 
a "counting device" . Indeed, the current-carrying wire may be replaced by a 
tube of curves, something that depends on v(r) alone: To this end, consider 
some closed curve C in the region where v(r) is defined. Consider further a 
family of integral curves of v(r) passing through C. These curves constitute a 
tube of curves. Now, the flux of v(r) across the surface S bounded by C remains 
constant along the tube. That is, /„ v • ridS — /„, v • n'dS' for all surfaces S' 
that cut the tube. This result follows from V • v = alone. Let us refer to the 
flux (J = J„ V • ndS* as the strength of the tube [13]. It is always possible - with 
an appropriate redefinition of v, if necessary - to let the tube of curves shrink 
into a line, keeping fixed the value of a. It is clear that a can play just the same 
role that the current / has played before. That is, we may substitute v(r)(i^r 
by crdx, with dx a length element of the curve x(s), and then proceed further 
as we did before: equation (13) reduces in this case to 

o" /" , (y ^ r) 
drx ^ ^ 

'ci |y-r|' 

so that the right sides of equations (16) and (17) acquire a factor a. If the 
curve Ci : x(s) crosses the surface 5* enclosed by C2 a number n of times, it 



W(y) = ^ / drx i^—^, (19) 



contributes with na to the total flux through S that appears on the left side 
of equation (17). The common factor a cancels then, and we arrive at Gauss' 
result again. However, this time the strength - at variance with the current - 
does refer to a geometrical property alone. 

By introducing the strength we are provided with a purely geometrical con- 
cept on which our proof of Gauss' result can rely. Nevertheless, for the sake of 
demonstrating Gauss' result, the introduction of the strength appears still as a 
disposable artifice. As already said, it is in fact possible to demonstrate that 
Jg V • ndS = n without introducing the strength as a substitute of the current, 
as it is done in the Appendix. 

In any case, we have arrived at Gauss' result by applying purely geometrical 
facts. We do not claim that the present method resembles more likely what 
Gauss probably did, than the method conjectured by the editor of Gauss' Works. 
Although it is true that the genius of Gauss could have anticipated some of 
Hclmholtz's results needed for the proof given above, there are other aspects 
that one should take into account before making any conjecture about the line 
of reasoning actually followed by Gauss. We are not concerned, however, with 
this kind of issues here. 

Our analysis has brought into light that the proof of Gauss' result that 
was based on Biot-Savart's and Ampere's laws was in fact based on Maxwell's 
equation V x (y^^^B) = J. We also saw that this last equation could be replaced, 
for the sake of proving Gauss' result, by the equation V x W = v. This equation 
can be seen as the one defining a vector field W(r), given another vector field 
v(r) whose divergence vanishes: V • v = 0. This one restriction put on v(r) is 
a particular instance of the continuity equation. We are so naturally led to ask 
about the equation that would replace V x W = v, when instead of V • v = 
we require from v to satisfy the continuity equation: dtp + V • pv =0. Putting 
it otherwise, we are led to ask about the connection between the continuity 
equation and Maxwell's equations. This is the issue we want to discuss in what 
follows, with the help of Helmholtz's theorem. 

4 Helmholtz's theorem 

In this Section we review the aforementioned Hclmholtz's theorem [8, 9, 10, 11, 
12]. According to this theorem, a vector field v(r) is completely determined by 
its divergence and its curl, together with a boundary condition which specifies 
the normal component of the field, n • v, at the boundary of the domain where 
the vector field is to be determined. For physical applications it is natural to 
take as "boundary" the infinity and the vector field vanishing there. Helmholtz's 
theorem then says that wc can write v(r) in terms of two potentials, ip and A, 
as 

v(r) = -V(p + V x A, (20) 

where (p and A can be expressed through the divergence and the curl of v(r). 
Writing 

10 



V • v(r) - p{v) (21) 

Vxv(r)=j(r) (22) 

we have that 

v(rO = -v/dV.-^^ + Vx ld\.-p}—^. (23) 

J 47r|ri-r2| J 47r |ri - rg 

This can be written, alternatively, as 

v(ri) = V I d\2G{Yi,Y2)p{v2) - V X /"dV2G(ri,r2)j(r2), (24) 

with the Green's function G'(ri,r2) satisfying 

V2G(ri,r2) = (5(ri-r2) (25) 

and vanishing at infinity. Under this last condition, the solution of equation 
(25) is given explicitly by 

Now, assume that we prescribe only the divergence p of a field, which is 
a function not only of position but of time as well, i.e., p ~ p{t,v). Let our 
boundary condition be that p vanish at infinity. Hclmholtz's theorem states 
that there is a field, call it D, such that 

V-D(t,r)=p(t,r). (27) 

The field D(t,r) is explicitly given by 

D(t, r) = -Vr / d^n-^^l^^ + Vr X Zi, (28) 

J 47r |r — ri 

with Zi an arbitrary field that we may take equal to zero, for simplicity. Note 
that the time t plays, in all of this, only the role of a parameter that can 
be appended to the fields, without having yet any dynamical meaning. The 
field D(t, r) satisfies therefore the only condition we have put upon it, i.e., 
V • D(t, r) = p{t, r). Its curl has been assumed as unspecified or else set equal to 
zero. We have then, by applying the gradient operator that appears in equation 
(28), 

D(i,r) = ^/rfV,^i^il^. (29) 

47r7 |r-ri|' |r-ri| 

For a point-like charge we set p(t, r) = q5{Y — YQ{t)) and the above expression 
reduces to 



11 



D(t,r) = ^ ,^-^. (30) 

According to equations (29) and (30), the field D(i, r) at time t corresponds 
to an instantaneous Coulomb field produced by a charge distribution p, or else 
by a point-like charge q. Such a result would correspond to an instantaneous 
response of the field to any change suffered by the charge distribution. Appar- 
ently, there is here a contradiction with the finite propagation-time needed by 
any signal. This issue has been discussed and cleared, in the case of the com- 
plete set of Maxwell's equations, by showing that both the longitudinal and the 
transverse parts of the electric field contain such instantaneous contributions, 
which turn out to cancel each other [14]. Note that by taking Zi equal to zero in 
equation (28) we have V x D = in our case, which is not what happens when 
D has to satisfy (together with H) the complete system of Maxwell's equations. 

Coming back to Hclmholtz's theorem, assume next that p(t, r) satisfies the 
continuity equation 

dtp + V-i^O. (31) 

By using equation (27) the continuity equation can be written as 

W-{dtB+j)=0. (32) 

Here, again, we can apply Hclmholtz's theorem and write the field dtT> + j 
as the curl of some other field H, whose divergence we do not need to specify. 
This field is given by Hclmholtz's theorem as 

H(^,r)=VrX fd^r,^^^hlllp^(^ + S7^, (33) 

J 47r|r-r2| 

As before, we may take V • H = 0, and hence ip = 0. That such H(i, r), with 
if set equal to zero, fulfills the required equation, namely 

V X H = j + dtB, (34) 

can be shown as follows. By writing 

A«,,./.=.M(ti|I±*^ ,35, 

in equation (33), we have H(t, r) = Vr x A(i, r), so that taking the curl on both 
sides of this last equation we obtain 

Vr X H(t, r) = V^ X Vr X A(t, r) = Vr(Vr • A) - V^A. (36) 

The first term on the right-hand side vanishes: 

Vr(Vr-A) = 0. (37) 
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This can be seen by applying the operator Vr(Vr-) to A given by equation 
(35). One replaces the second derivatives with respect to r in the term |r — r2|~ 
inside the integral, by second derivatives with respect to r2. Integrating by 
parts one obtains equation (37), upon using that the field j(i, r) + 9fD(i,r) 
is solenoidal, i.e., one whose divergence vanishes, and that it is bounded in 
space, or that it vanishes faster than r~^ for large r [12]. We have therefore 
V X H = -V^A, so that 

V.xH = -Jd\2m,r2)+dtB{t,r2))Vl(—j^—^Y (38) 

Taking now into account that 

1 



wc finally obtain 



V^(t^-— )=-4^^(r-r2), (39) 



VxH=j(f,r)+atD(t,r), (40) 



which is identical to one of the two Maxwell's equations with sources. The other 
Maxwell equation that has a source term is identical to equation (27). 

5 Maxwell's equations with sources 

Let us summarize how the inhomogeneous Maxwell's equations followed from 
the assumption of charge conservation. Let us describe charge - or any other 
quantity, e.g., matter - by its density p{t, x) and its velocity distribution v(i, x). 
It is well known that charge or matter conservation can be expressed through 
the equation of continuity 

dtp + V- (pv) = 0. (41) 

As we have shown above, given a bounded, scalar, function p(i,x), there 
exists a field Y(i,x), satisfying 

V ■ Y = p, (42) 

so that equation (41) can be written in the form 

V • {dtY + pv) = 0. (43) 

Because the field dtY+pv is divergenceless, there exists another field W(i, x) 
such that 

V X W =dtY + pw. (44) 

Equations (42) and (44), written with the usual notation (Y ^ D, W ^ H, 
pv^ j), i.e., 
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V • D = p (45) 

V X H = j + dtT>, (46) 

are the two inhomogeneous Maxwell's equations. They follow from equation 
(41) alone. Reciprocally, the continuity equation follows from equations (45) 
and (46) by taking the divergence of the second of these equations and replacing 
in it the first equation. We stress that the curl of D and the divergence of H 
have been left unspecified - or set arbitrarily equal to zero - when going from 
equation (41) to equations (45) and (46). In the case of the electromagnetic 
field, such an arbitrariness does not correspond to the observed facts, as charge 
and current distributions completely determine the fields they produce. There 
must be, therefore, additional equations for the curl of D and the divergence of 
H in order to completely determine these fields, when their boundary conditions 
are also specified. 

Let us now turn our attention to the source-free, or homogeneous. Maxwell's 
equations: 



V ■ B = (47) 

V X E+9tB = 0. (48) 

Of course, by a similar reasoning as above, we could derive these equations 
from the conservation of something, call it pm- The conservation of pm could be 
expressed through an equation similar to the equation of continuity, equation 
(41). After deriving equations similar to (45) and (46) we put pm — and 
obtain equations (47) and (48). In other words, this something - a "magnetic 
charge" if you want - is assumed to be conserved: zero before and zero after. 

Alternatively, the source-free Maxwell's equations can be understood as a 
purely mathematical statement telling us that the fields E and B derive from 
a scalar function ip and a vector potential A (not to be confused with those in 
equation (20)), through the definitions 



E = -Vip - dtA, (49) 

B = V X A. (50) 

In other words, given two arbitrary fields, (p and A, by defining E and B 
through equations (49) and (50), it follows that equations (47) and (48) are 
identically satisfied. 

In any case, we should stress that equations (45, 46) on the one hand, and 
equations (47, 48) on the other hand, are - up to this point - independent from 
each other. We may connect them through some constitutive equations, like, 
e-g-, 
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D = eE (51) 

H = ^"^B. (52) 

These equations are usually assumed to describe a linear medium of electrical 
permittivity e and magnetic permeability /i. A particular case of such a medium 
is vacuum, and the system of equations (45, 46, 47, 48), that arises out of a 
connection like the one given by equations (51, 52) is what we know as the 
complete system of Maxwell's equations. 

Let us stress once again that without connecting (D, H) with (E, B) through 
some constitutive equations we have no closed system. The equations that we 
have written down for (D, H), that is Maxwell's equations with sources, can 
also be written down for a fluid, for example. These equations are, as already 
stated, V- Y = p and V x W = j+i9tY. We can expect that any conclusion that 
can be derived in the realm of electrodynamics from the equations V • D = p and 
V X H = j + dtT) without coupling them to the source-free Maxwell's equations, 
will have a corresponding result in the realm of fluid dynamics. This assertion 
can be illustrated by two examples: 1) The case in which a point-like singularity 
of p is produced in the interior of a fluid, for which one obtains a velocity- field 
obeying a law that is mathematically identical to Coulomb's law [13]. 2) The 
case of a fluid where a so-called vortex tube appears (tornadoes and whirl-pools 
are associated phenomena) , in which case - after approximating the vortex-tube 
by a line singularity - one obtains a velocity-field through an expression which 
is mathematically identical to the Biot-Savart law [13]. 

Finally, let us mention that the derivation of the inhomogeneous Maxwell's 
equations as a consequence of charge conservation is nothing new [15]. It follows 
as a direct application of a theorem of de Rahm for differential forms: given a 
4- vector j°'{x) for which a continuity equation holds, 9a j" = 0, there exists an 
antisymmetric tensor F""^ — —F^'^ such that daF"'^ = j^. This last equation 
is nothing but the tensorial form of the inhomogeneous Maxwell's equations 
(45) and (46). Now, the tensor F"^ does not need to be derivable from a 4- 
vector A", according to F^^ = d°'A''^ - S'^A". For this to be the case, F"!^ 
would have to satisfy d°'F'^'^ + d'^F'^'^ + ^tF"'^ == 0. This is the tensorial 
form of the homogeneous Maxwell's equations (47) and (48). In other words, 
given j" and ^", with j" satisfying a continuity equation, we can introduce two 
antisymmetric tensors, F?"^ and F'^^, the first one satisfying daF?^ = j^ and 

the second one, F^'^J = 9"A^ - d'^A°', satisfying ^"Fg]' + d^F^^^ -h d-<F^^^ = 0. 
In order that these two equations do conform a closed system, i.e. the total 
system of Maxwell's equations, we need to connect ^"f with F°Z through some 

constitutive relation®. 



®As a matter of fact, it is not necessary to rest on de Rham's theorem and the theory 
of differential forms on manifolds in order to derive the foregoing conclusions in tensorial 
form. Indeed, one can start with the tensorial form of Helmholtz's theorem (see Ref. [9]) and 
go-ahead with a similar reasoning as the one we have followed in the preceding paragraphs. 
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6 Suminary and Conclusions 

We have presented an elementary proof of Gauss' result concerning the linking 
number of two curves. The proof is based on geometrical properties of a vector 
field, properties that are at the root of Helmholtz's theorem. It could be cleared 
up what at first glance appeared as an astonishing connection between Gauss' 
linking integral of two curves and the line integral of a magnetic field. By ex- 
posing the connection that holds between Helmholtz's theorem and the linking 
integral of Gauss, we could show which purely geometrical properties were be- 
hind the line integral of a magnetic field. In this way, we could lay bare the 
essential constituents of the connection between an empirically established nat- 
ural law (Maxwell's equations) and a geometrical property of two curves (their 
linking integral). Such essential constituents are certainly not characteristic of 
the electromagnetic field. As a consequence, we showed that starting from the 
conservation of any property - like mass or charge - for which the continuity 
equation holds true, one can arrive at equations which are mathematically iden- 
tical to the Maxwell's equations. These equations are therefore tightly linked to 
a general statement telling us that something is conserved. Consider anything 
- charge, matter, or whatever - the amount of which that is contained inside 
an arbitrary volume changes with time. If this change is exclusively due to a 
flow through the volume's boundary, then a continuity equation holds true. As 
a consequence of it, a pair of Maxwell-like equations must be fulfilled for some 
auxiliary fields that we may introduce, playing the role of D and H. 

The fact that Maxwell's equations with sources can be derived from charge 
conservation is known [15], although it is not usual to find it mentioned in stan- 
dard textbooks of electromagnetism. Maxwell's equations with sources involve 
the fields D and H, whereas the source-free equations involve the fields E and 
B. It is through some constitutive equations connecting (D, H) with (E, B) 
that we obtain a closed system, i.e., the complete system of Maxwell's equa- 
tions. The constitutive equations express, in some way, the underlying proper- 
ties of the medium where the fields act or are produced. From this perspective. 
Maxwell's equations entail besides charge conservation some other properties of 
the medium, yet to be unraveled. These properties are effectively described, in 
the simplest case, through the permittivity e and the permeability /i. The first 
one refers to electrical, the second one to magnetic, properties of the medium, 
be it vacuum or any other one. It is just when the equations for (D, H) to- 
gether with those for (E, B) do form a closed system, that we can derive a wave 
equation for these fields. The velocity of wave propagation is then given by 
c = (e/i)^^/^, the velocity of light. It is remarkable that the velocity of light can 
be decomposed in terms of a product of two independent parameters. However, 
the historical development in physics has led us to look at c as a fundamental 
constant of Nature, instead of e and /i . Nevertheless, currently discussed and 
open questions related to accelerated observers, Unruh radiation, self-force on a 
charge, magnetic monopoles and the like, might well require an approach where 
the role of c recedes in favor of quantities like e and fi. Such an approach could 
perhaps be reminiscent of Helmholtz's efforts to develop an ether theory as the 
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basis of electrodynamics, a theory that strongly relied on the mathematical ma- 
chinery that Gauss so much helped to develop. If this happens to be the case, it 
will seem somewhat ironical that the Maxwell's equations, when written in the 
- by now - most commonly used Gaussian units, do not include but the single 
constant c, hiding so e and /x from our view. These last two constants might 
well be key pieces that remain buried under the beauty of a unified theory of 
electromagnetic phenomena, which is the version of electrodynamics that we 
know and use today. 

7 Appendix 

The assertion that Jg v • ndS = n can be proved in the following way, in which 
we will try not to veil with abstract definitions and too much mathematical 
detail what wc do intuitively take for true. 

The curve Ci pierces the surface S at different points Pi , P2, ■ ■ ■ , Pn ■ These 
points can be obtained as intersections of Ci and S"^ . Consider a neighborhood 
Mi of Pi, it being sufficiently small so as to contain no other point Pj ^ Pi. We 
can deform this small portion A/i of S, making it flat and lying parallel to the 
XF-plane. By so doing we do not change the value of /^ v • ndS, as long as C2, 
the bounding curve of 5, remains unchanged. The idea is to split the integral 
/„ V • ndS into N contributions, each one of them corresponding to a small, flat 
surface that locally can be made lye parallel to the XY-plane. Now, in general, 
the integral Jg v • ndS can be expressed as [7] 



v-nd5= / Vz{x,y,f3(x,y))dxdy. (A.I) 

s Js 

Here, v = {vx,Vy,Vz), whereas S is assumed to be described by the equation 
z = fi{x,y). With V given by equation (14) and after properly adapting the 
variables in equation (A.I) to our case, we have 



V ■ ndS = ds—S{rx - x{s))S(ry - y{s))S(rz - z{s))drxdry. (A. 2) 

In order to evaluate the integral on the right-hand side of equation (A. 2) we 
use the rule 

"^ \9{xi)\ 

where the Xi satisfy g{xi) — 0. Applying this property successively to equation 
(A. 2) we obtain 

^Let Ci be given by equations of the form x' = /*(s), and 5 by equations of the form 
x^ = g^{u,v) with i = 1,2,3 and s,u,v some parameters. The points of intersection Pj. G 
Ci n S are obtained by solving the system f'{s) = g^{u, v), i = 1,2, 3, for s, u, v. In general, 
there are a number Af > of solutions (s],,U],, vi^), k = 1, . . . N, oi these equations. 
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/ v-ndS^'Y] / jda—5{r,^-x{s))5(ry-y(s))5(r^_(rx,ry)-z(s))dr^dry 
Js i^i J^f^ J "* 

N „ N 

= ^ ± / 5{rx - x{s))S(ry - y(s))dra:dry = ^ ±1 = n < TV. 

Z— 1 -^^ 2—1 

We have used the fact that the surface integral sphts into N integrals, i.e., 
as many as the number of times that Ci pierces 5. We have also used that each 
Afi can be described as a small part of the plane Vz — const.. It might occur 
that in some cases Ci goes through 5* in one direction, turns then back and goes 
through 5* again but in the opposite direction, without having wound around the 
boundary C'2 of S in-between. For these cases we have a positive contribution 
(dz/ds — + \dz/ds\) followed by a negative one (dz/ds — — \dz/ds\). These 
contributions cancel in pairs. In order that the contributions of two consecutive 
points Pi and P^+i are of the same (positive) sign, it must occur that the portion 
of C'l joining Pi and P^+i winds once around C2. At the end, what we obtain as 
a result of the sum X]i=i(^l) i^ ^^^^ number n of times that Ci winds around 
C2. Hence, we may use that Jg v • lids' = n in equations (17) and (18) to obtain 
Gauss' result. 
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